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5 METHOD AND APPARATUS FOR IMPROVED DEPTH 

MATCHING OF BOREHOLE IMAGES OR CORE IMAGES 

Field of the Invention 

10 This invention relates to logging of boreholes. More specifically, the invention relates 

to a method for improving correlation of core images, borehole images and well logging data 
derived from sensors located on a well logging tool as it passes through a borehole. 

Background of the Invention 

15 

In a typical prior art oil well logging scenario, a string of well logging tools having 
multiple sensors for measuring characteristics of the earth formation along the wall of a 
borehole is lowered via a cable to the bottom of the borehole. Geophysical data is recorded 
by way of the sensors as the cable is wound in using a precision winch. The depth in the 
20 earth at which the sensors on the logging tools are positioned as data is logged is determined 
by measuring the logging speed and cable depth. Devices such as a depth wheel 
measurement instrument and an axial accelerometer may be utilized. 

Typically, measurements taken along the length of a borehole by logging tools are 
25 intended to provide indications of oil-bearing or gas-bearing strata in the earth. In the prior 
art, measurements of various characteristics or parameters of earth formations are usually 
obtained by combining measurements (data logs) taken by way of multiple sensors during a 
single pass through the borehole, or taken during different passes through the same borehole. 
When combining such measurements it is necessary that they be accurately correlated in 
30 depth with one another to be useful. 

One problem is that the movement downhole of a tool with its sensors is usually not 
uniform. The non-uniform motion may be caused by such thing as: (a) damped longitudinal 
oscillations of the logging tool on the cable, (b) sticking and slipping of the logging tool 
35 against the sides of a borehole or the wall of casing in the borehole, and (c) irregular motion 
of the sensors that are mounted on mechanical arms that have independent motion with 
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5 respect to the logging tool itself. As a result of non-uniform motion of the logging tool, data 
collected by a particular sensor on the logging tool at any specific depth in the borehole may 
not show as being recorded at the specific depth but at a different depth. Similarly, data from 
different sensors of the same tool may not show as having been recorded at the same depth 
for all those sensors. 

10 

To overcome the depth measurement problem in the prior art, a well logging tool for 
measuring depth in a borehole is disclosed in U.S. Patent 5,019,978, issued May 28, 1991 to 
Allen Q. Howard, Jr. and David J. Rossi. The invention in that patent provides for estimation 
of a dominant mechanical resonant frequency parameter and of a damping constant 
15 parameter. These two parameters are taken into consideration when correcting an 

approximate indication of depth of a well logging tool to determine the actual, true depth of 
the well logging tool in a borehole. 

Besides measuring depth in a borehole, well logging seeks to measure other physical 

20 characteristics along the sides of the borehole such as fractures, bed boundaries and bed dips. 
A major advance in borehole logging has been the development by Schlumberger of the 
Formation Microscanner ("FMS"), a borehole imaging system. The FMS processing 
technique is described in U.S. Patent Numbers 4,468,623 issued Aug. 28, 1984 to Stanley C. 
Gianzero, David E. Palaith, and David S. K. Chan; and in U.S. Pat. No. 4,567,759 issued Feb. 

25 5, 1986 to Michael P. Ekstrom and David S. K. Chan. The FMS system uses a tool having 
wall-engaging pads each carrying an array of electrical sensors distributed in the 
circumferential direction with respect to the axis of the borehole. Signal voltages generated 
by the sensors are sampled as the well logging tool moves along the borehole. The signals 
are processed and rendered visible, by photographic or other printout, or by cathode-ray tube 

30 display as a two-dimensional visible image is formed over logged segments of the borehole 
walls. In such images, bed boundaries can be visually identified from sharp visible contrasts 
in the images, which reflect sharp changes in resistivity at boundaries of the beds. The 
images thus obtained may exhibit a resolution on the order of 0.5 cm, allowing very fine 
details of the formation to be distinguished due to the number of sensors in the 

35 circumferential direction, and the high rate of sampling in the longitudinal direction. 
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5 U.S. Patent No. 4,251,773, issued Feb. 17, 1981 to Michel Cailliau and Philippe 

Vincent teaches the use of signals from sensors to determine dip (inclination) and azimuth 
(strike) of bed boundaries. More specifically, this patent teaches a logging tool that has four 
substantially identical pads with sensors angularly distributed about the axis of the logging 
tool in a side-by-side relationship and adapted to engage the borehole wall at ninety-degree 

10 intervals. The sensors provide resistivity measurements of the respective sectors of the 

borehole wall engaged by the pads. As the logging tool is moved along the borehole wall, the 
sensors continuously provide signals measuring the resistivity of the adjacent earth formation. 
Sharp variations in resistivity indicate boundaries between different beds in the earth 
formation. The signals produced by the sensors at different angular positions of the pads are 

15 processed to provide information about the dip of bed boundaries, i.e., the orientations of the 
bed boundaries with respect to a terrestrial reference, and the azimuth of the dip. 

Another tool that is lowered into or withdrawn from a borehole and has multiple 
sensors that measure properties such as resistivity along segments of a borehole wall is taught 

20 in U.S. Patent 5,960,371 issued September 28, 1999 to Naoki Saito, Nicholas N. Bennett and 
Robert Burridge. This patent teaches use of the Hough transform to extract dip and azimuth 
of multiple fractures and beddings from any type of borehole image with respect to a 
terrestrial reference. The method is also robust enough to account for noise or gaps in the 
images. The method can separate dips and azimuths of fractures from those of formations. 

25 Thus, it can detect and characterize other geometric features (e.g., linear, circular, or 

ellipsoidal shapes, some of which may represent vugs in carbonate reservoirs) present in the 
images. 

As previously pointed out, when combining data of well logs it is necessary that they 
30 be correlated in depth in order to be useful. One method for depth correlation of well log 
data is taught in U.S. Pat. No. 4,327,412 issued April 27, 1982 to John P. Timmons. The 
method disclosed in this patent is relatively complex and determines the displacement 
between a plurality of well logs so they may be correlated and combined. The well logs are 
derived from multiple, spaced sensors passed one time through a single borehole, or from 
35 separate passes of the same sensors through the same borehole. First, a normalized 

correlation function between selected groups of samples of the sets of data is determined as a 
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5 first assumption to have a predetermined displacement relationship with the groups of 

samples of the data. A step of determining the normalized correlation function is repeated for 
a number of overlapping groups of well log data samples to produce a number of overlapping 
correlation functions. At least some of these overlapping correlation functions are combined 
to produce an improved correlation function that is used to depth correlate the well log data to 
10 a common, accurate depth level. 

Another method for correlating well logging data collected by multiple sensors is 
taught in U.S. Patent 6,272,232 issued August 7, 2001 to Jean-Pierre Delhomme and Jean F. 
Rivest. This patent teaches a method for constructing, from an initial image of the wall of a 

15 borehole, a new "crossing-component image" centered on the axis of the borehole. The new 
image is representative of variations in a physical parameter of the earth formation in both the 
longitudinal direction of the borehole (depth), and in the peripheral direction of the borehole 
wall (laterally). The new image includes only those components of the physical parameters 
that extend all the way across the initial borehole image from one side of the image to the 

20 other. The method also includes determining variations in one or more attributes relating to 
the new image as a function of depth. The variations provide information relating to 
morphology to indicate solid zones, bedded zones, or different types of heterogeneous zones. 

When individual well logs are concurrently obtained from a first and a second set of 
25 sensors vertically spaced from each other in a well logging tool the well logs also need to be 
depth correlated. They need to be depth correlated so that the data for any given level in the 
well taken by the two sets of vertically spaced sensors are aligned in order to be useful. 
Another method for depth correlating well log data is taught in U.S. Patent 4,320,469, issued 
March 16, 1982 to William J. Frawley and Philip A. Mongelluzzo. This patent teaches doing 
30 this by using data correlograms obtained by applying a correlation function to a pair of 

digitized well logs to reduce the amount of data that must be processed. The remaining data 
is then processed in an efficient and accurate manner in order to arrive at results indicating 
exactly how much two logs must be shifted with respect to each other for an optimized depth 
correlation between them. 
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5 While methods for depth correlating data obtained using sensors spaced vertically 

from each other on a well logging tool, or obtained at separate times of the same borehole are 
known, some of these methods are very complex and there is a need for a better method to 
depth correlate data obtained from well logging tool sensors. 

10 Description of the Drawings 

The invention will be better understood upon reading the following Detailed 
Description in conjunction with the drawings in which: 

15 Figure 1 shows a well logging tool in a borehole along with apparatus at the surface of 

the earth for controlling the tool and recording measurements derived therefrom; 

Figure 1 A is a highly simplified representation of a portion of a well logging tool 
showing the orientation of sensors thereon; 

.20 

Figure 2 is a picture of a portion of a well logging tool with its arms, pads and 
sensors; 

Figure 3 A shows a bedding plane or fault intersecting a borehole at an angle relative 
25 to the axis of the borehole to create a dipping event and the intersection is represented by an 
ellipse; 

Figure 3B shows a representation showing a portion of the wall of a borehole 
represented in Figure 3A unrolled to a flat position; 

30 

Figure 4 shows a display of signals obtained from sensors mounted on two sets of 
pads angularly distributed about the axis of a well logging tool, with each set of pads being 
spaced vertically with relation to each other; 

35 Figure 5 shows signals obtained from two sets of sensors spaced vertically with 

relation to each other after depth correlating the signals for vertical offset; 
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Figure 6 shows a block diagram of the steps involved to perform depth correlation by 
averaging individual sensor signals obtained at each level of a well logging tool into a single 
signal, and the averaged signals are then depth correlated to determine the optimum depth 
offset necessary for precisely matching signals; 

10 

Figure 7 shows averaged two dimensional signals obtained from two sets of sensors 
spaced vertically with relation to each other after being translated into three dimensional 
Hough space to be used to determine depth correlation to be applied to the signals; and 

15 Figure 8 shows a block diagram of the steps involved in translating signals taken from 

the two sets of sensors into three dimensional Hough space to perform depth correlation. 

Detailed Description 

20 The terms "depth match" and "depth matching" used hereinafter mean correlating sets 

of data taken by different vertically spaced sensors belonging to different levels of sensors, 
such as sensors S la - Sid and S2a - S2d in Figure 1A, during a single pass of a logging tool 
through a borehole, or during different passes through a borehole. For a geological event, the 
vertical offset between different levels of sensors causes signals coming from sensors in a 

25 first level of a logging tool to be generated at a slightly different time than the signals coming 
from sensors in a second level. The difference in these times is also dependent on logging 
speed variations. The effects of the time-offset need to be compensated to "depth match" the 
signals coming from the sensors at both levels. 

30 Briefly, in a first embodiment of the present invention, individual traces from the 

sensors of any first specific level of sensors on a logging tool are converted into one signal 
representing those traces. The one signal represents all the traces (averaged log signal) and 
includes an average borehole-crossing signal. The averaged log signal is generated by 
performing a prior art computation of bedding dips or flowlines (angle and azimuth) on the 

35 traces from each sensor and an average or median computation is applied to all the samples 
from sub-sensors that comprise each sensor aligned along the pre-determined dip. This is 
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5 repeated for a second specific level of sensors vertically offset from the first specific level of 
sensors and a second averaged signal is obtained. Then a prior art single log trace depth 
matching technique is applied to two averaged signals to determine the optimum depth 
"offset" necessary to precisely match the two averaged signals. The "offset" determined is 
then applied to the original traces from the sensors at one level to depth match them to the 
10 traces from the sensors at the other level. 

Briefly, in a second embodiment of the present invention the Hough transform is 
utilized to generate three-dimensional images from well log signals. The three dimensional 
images are used to create two-dimensional extremum curves that are in turn used to calculate 

15 an offset. The offset is used to depth match the well log signals obtained from vertically 
spaced sensors on a well logging tool, or signals obtained during different logging passes 
through the same borehole. The detrimental effects of partial and non-overlapping sensor 
coverage of a borehole wall are reduced, if not eliminated, because of the characteristics of 
the Hough image transform, and because depth matching is finally applied to the two 

20 dimensional extremum curves derived from the three-dimensional images. More specifically, 
a depth matching offset is calculated directly from two dimensional extremum curves derived 
from the Hough three-dimensional images without the need for calculating dip and azimuth 
of bed boundaries in an earth formation and compensating for it before the depth matching 
may be performed. 

25 

While the following description is given with reference to depth matching data 
obtained from sensors on a borehole logging tool traveling through a borehole, the teaching 
of the present invention may be utilized with records, no matter how obtained or made. The 
records may be a large number of azimuthally oriented traces obtained either from sensors 
30 applied to a borehole wall, from sensors scanning a borehole azimuthally, from sensors 
scanning azimuthally the cylindrical surface of a whole core, or from sensors probing an 
azimuthally oriented core section cut, by stacking the pixels/measures along the apparent dip 
directions seen on the core section cut. 

35 In the following description only two sets of data are obtained from two sets of spaced 

sensors, and the data is correlated. This is done for the sake of clarity in describing the 
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5 invention. However, the teaching of the invention may be used to correlate or depth match 
multiple sets of data. In addition, while the sets of data collected are described as being two- 
dimensional, some types of sensors may collect three-dimensional data. 

In the following description of a well logging system, prior art material is not given in 
10 great detail to avoid detracting from an understanding of the two embodiments of the 

invention presented herein. For a more detailed understanding of any referenced prior art see 
the above Background of the Invention and the references cited therein. 

Referring to Figure 1, there is shown an investigating tool 10 in a borehole 11 for 
15 investigating subsurface earth formations. Investigating tool 10 is supported in borehole 11 
on an end of a cable 12. Cable 12 passes over a sheave wheel 13 and is secured to a drum 
and winch mechanism 14. The drum and winch mechanism 14 includes a suitable brush and 
slip ring arrangement for providing electrical connections between the cable conductors and a 
control panel 15. The control panel 15 operates to supply power and control signals to 
20 investigating tool 10, and includes suitable electronic circuitry for receiving well logging 
signals from sensors on the investigating tool. The well logging signals may be stored on a 
digital tape recorder 16. 

Tape recorder 16 converts analog signals received from control panel 15 into digital 
25 signals and is stepped as a function of depth by a driving wheel 17 that engages cable 12. 
The digital signals are then either transmitted or carried to a computer system 18 for 
processing. Computer system 18 comprises a processor, display and plotter for processing 
and displaying well logging data. 

30 Investigating tool 10 has a reference point 20 that constitutes the center or recording 

point of tool 10. Because of long elastic cable 12, investigating tool 10 is subject to 
displacement arising from cable stretch thus causing the true depth level to be different from 
the depth level registered by a depth recorder 19 derived from driving wheel 17. The depth 
level registered by depth recorder 19, which is identical with that registered by the tape 

35 recorder 16, is designated "Z". When attempting to depth match well log data produced by 
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5 different sensors of investigating tool 10 there is a distinct possibility that the well logs 
recorded at the same depth level will not be referenced to the same depth level. 

On investigating tool 10 sensors are located on two vertically spaced levels LI and L2 
as shown, and their positions are described in greater detail starting in the next paragraph. 
10 These sensors produce the well logging signals that are received by control panel 15. 

Figure 1A is a simplified representation of a portion of tool 10 about reference point 
20 showing the layout of sensors thereon. Sensors are usually mounted on pads that engage 
the wall of a borehole or casing in a borehole, although the pads are not shown in Figure 1 A. 
15 In the highly simplified arrangement described herein there is only one sensor on each pad. 
In Figure 1A the sensors are represented by boxes with an "x" therein and are designated Sla 
- Sid in level LI and S2a - S2d in level L2. In Figure 2 there is shown a more detailed 
representation of an investigating tool 10. 

20 Figure 1 A shows two levels LI and L2 of sensors positioned in two parallel planes. 

Level LI has four sensors Sla - Sid spaced ninety degrees apart, and level L2 has four 
sensors S2a - S2d spaced ninety degrees apart. In one embodiment, the level L2 sensors are 
rotated or offset forty-five degrees with respect to the level LI sensors. The spacing of the 
sensors in each of levels LI and L2, and the rotational offset between the two planes, results 

25 in the orientation of all the sensors around a borehole as represented on the circle at the top of 
Figure 1A. 



An example of a sensor that may be used on tool 10 is one that measures resistivity or 
conductivity to detect and measure boundary layers in a borehole. Such a sensor is taught in 
30 U.S. Patent 4,786,874 issued November 22, 1988 to Donald S. Grosso and Allen Duckworth. 
The patented device utilizes electrical conductivity or resistivity to detect and measure 
boundary layers. 

Figure 2 is a picture of an actual well logging tool 10 with its arms 52a - 52d, pads 
35 53a - 53d, and sensors Sla - Sid thereon. Arms 52a - 52d move independently as tool 10 
moves through a borehole to account for irregularities in the wall of the borehole. These 
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5 arms also move perpendicularly to the axis of a borehole to keep sensors on the pads in 

contact with the wall of the borehole. Pads 53a - 53d are respectively connected to arms 52a 
- 52d and move therewith. On the surface of pads 53a - 53d are mounted sensors Sla - Sid. 
Each of the sensors Sla - Sid may include a number of smaller sub-sensors. 

10 When a bedding plane or fault intersects a borehole 1 1 at an angle relative to the axis 

of the borehole, it is inclined or dipping and creates a "dipping event" and the intersection is 
represented by an ellipse 57 as shown in Fig. 3a. In Fig. 3A are representatively shown four 
traces 2 1 a - 2 1 d taken by the four sensors Sla-Sldin level L 1 of Fig. 1 A as the tool moves 
through the borehole. The traces are not shown extending longitudinally along the periphery 

15 of borehole 1 1 to avoid cluttering up Fig. 3 A. The four traces intersect ellipse 57 at four 

points Zi - Z4. These four points represent the different times at which the four sensors Sla — 
Sid will generate a signal representing the same dipping event as the sensors in level LI are 
passed through borehole 11. 

20 When the portion of the wall of borehole 1 1 represented in Figure 3 A is unrolled to a 

flat position, the ellipse 57 appears as a sinusoid 58 as shown in Fig. 3B. The location of 
trace- intersect points Zi - Z4 are shown on sinusoid 58. In addition, if traces 22a - 22d in 
level L2 of Fig. IB were shown plotted on Fig. 3 A, which they are not since they are not in 
level LI, they would appear in Fig. 3B as the four intersect points Z 5 - Z 8 and are interspaced 

25 between points Zi - Z4 on sinusoid 58 as shown. With all eight intersect points shown in Fig. 
3B it can be seen how a dipping event is represented as a sinusoid 58. 

When the images of a dipping event are projected onto a flat surface, as shown in Fig. 
3B, the vertical displacement of points Zj (where i = 1 through 8) with respect to a point Zo 
30 will satisfy the mathematical relationship: 
8Z0 = r (cp,Zo + 8Z0) tan (a) sin (q> -p) 
where: 

Zo= depth of the center of the event along the longitudinal axis of the borehole; 
cp = angle between north and the imaged point; 
35 5Zo = Z = vertical distance between the center of the pattern at level Zo and a point Z; 
a,p = dip and strike angle of the dipping event; and 
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5 r(<p, Zo+5Zo) = radius of the borehole at point (cp,Zo + 5Zo). 

Thus, the vertical displacement distance 5Zo from point Zo to each of points Zi - Z4 in 
level LI may be calculated. 

10 r(cp, Zo+5Zo) for round cylindrical boreholes is constant and the borehole intersection 

is an ellipse. When the ellipse is unrolled, as described above, the pattern seen is a sinusoid. 
The amplitude of this sinusoid is a function of the dip angle, and its phase is a function of the 
orientation of the surface with respect to north (strike angle or azimuth). For a non- 
cylindrical borehole, the intersection is not an ellipse and the pattern is a distorted sinusoid. 

15 Obviously, the amount of distortion depends upon the degree of departure of the borehole 
from a perfect cylinder. 

Fig. 4 shows displays 21 and 22 with traces 21a - 2 Id and 22a - 22d respectively. 
The set of traces 21a-21dis from the four sensors Sla-Sldin level LI, and the set of 
20 traces 22a - 22d is from sensors S2a - S2d in level L2 as seen in Figure 1 A. The traces are 
obtained during logging of borehole 1 1 . Each sensor may include at least one row of smaller 
sub-sensors. In addition, each sub-sensor generates a signal, and the side-by-side 
combination of such signals from the sub-sensors generates the eight traces 21a - 21d and 22a 
-22d. 

25 

It can be seen in display 21 that sensor Sla produces trace 21a at an arbitrarily 
assigned azimuth of 0 degrees. With the previously described ninety degrees spacing 
between sensors, sensor Sib produces trace 21b at an arbitrary azimuth of 90 degrees, sensor 
Sic produces trace 21c at an arbitrary azimuth of 180 degrees, and sensor Sid produces trace 

30 2 Id at an arbitrary azimuth of 270 degrees. The four level L2 sensors S2a - S2d that produce 
traces 22a - 22d have the same azimuthal spacing with respect to one another as the level LI 
sensors, i.e., the level L2 sensors are spaced apart at ninety degrees among themselves. 
However, as shown in Fig. 1 A, these sensors are oriented at a forty- five degree offset relative 
to the level LI sensors as explained with reference to Figure 1 A. Accordingly, trace 22a is at 

35 forty-five degrees, trace 22b is at one-hundred thirty five degrees and so on. 
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5 As logging tool 10 moves along borehole 1 1, it may rotate horizontally a small 

amount and its azimuth orientation may change. This creates the horizontal image skewing 
of the four traces shown in each of displays 21 and 22. The skewing is shown as the traces 
not being in a vertical, straight line. 

10 Superimposed on display 21 are drawn two sinusoidal waveforms 26 and 27 that 

reflect common borehole crossing features shown in traces 21a - 21d. "Borehole crossing 
features" designate rock features that are going across the borehole as opposed to smaller 
localized rock features that may be seen by fewer than all the pads of a level of sensors. 
Similarly, superimposed on display 22 are two sinusoidal waveforms 28 and 29 that reflect 

15 common borehole crossing features shown in traces 22a - 22d. 

As is known in the prior art, and as previously described with reference to Figs. 3A 
and 3B, sinusoidal waveforms 26 - 29 always reflect a "dip" in a bedding plane or fault that 
intersects borehole 11. The amplitude of the waveforms indicates the degree or angle of dip, 
20 and the location of the peaks of the waveforms (phase) indicates the direction or azimuth of 
the dip. 

As traces 21a - 2 Id and 22a - 22d of displays 21 and 22 are generated from sensor 
signals taken in the same borehole, they generally would look alike except for a vertical 
25 displacement 30, reflecting the spacing between the level LI and L2 sensors, and the azimuth 
offset between the level LI and L2 sensors. The vertical displacement of the traces in images 
21 and 22 reflects the fact that sensors S la - Sid are physically located above sensors S2a - 
S2d on logging tool 10 as shown on Fig. 1 A. In addition, some vertical displacement may be 
caused by other factors such as stretching of cable 12 as previously described. 

30 

Fig. 4 also shows single borehole crossing traces 23 and 24 that are obtained from the 
individual traces 21a-d and 22 a-d shown in this Figure. Each single trace 23 and 24 
respectively represents one set of traces 21a - 2 Id and 22a - 22d, each set being converted 
into single traces 23 and 24 respectively. In the first embodiment of the invention, traces 23 
35 and 24 are obtained by performing an average or median operation among each 
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corresponding set of traces 21a - 21d and 22a - 22d, each trace 23 and 24 representing the 
average value of a specific earth (rock) feature crossing the borehole. 



Each of the average traces 23 and 24 is generated by first performing a prior art 
computation of bedding dips or flowlines on the individual traces from each sensor of a 

10 specific level. As a result, sinusoidal waveforms, such as 26 are obtained. Then, an average 
or median computation is applied to all the sample images (responses from sub-sensors) 
located along the sinusoid waveform 26 (bedding dip) obtained. By averaging the signals 
along bedding dips, the common part of the signals that is due to bedding is enhanced and the 
localized part that is due to heterogeneities, such as gravel, vugs, etc., is attenuated. This is 

15 desirable as the common part may be expected to be found on all levels of sensors and thus 
be useable for depth matching correlation. 

The determination of a sinusoidal waveform relies on the dip and azimuth of each bed 
boundary. The dip and azimuth of a bed boundary as a function of depth in a borehole may 
20 be determined from traces of the borehole by using prior art techniques such as the Hough 
transform described in the background of the invention. 

As signals 23 and 24 are obtained, a single log trace depth matching technique is 
applied to these signals to determine an optimum depth "offset" for depth matching signal 24 
25 to signal 23. The "offset" determined is then applied to traces 22a - 22d to depth match them 
to traces 21a-21das seen in Fig. 5. 

The averaged traces 23 and 24 may be determined mathematically as follows. First, 
the intersection between a bedding plane and the borehole at a given depth Zo may be 
30 expressed mathematically as follows. 

Z = Zo + A(Zo) cos ( (p - (po(Zo)) 
where the Greek symbols cp and Z are the coordinate system of the borehole image space, and 
A(Zq) and (po(Zo) are the amplitude and azimuth of the bedding plane at the depth Zq. 
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5 The average signal trace at depth Zo, S(Zo), is then calculated as follows. Let I((p, Z) 

be the pixel values of the traces comprising a borehole image. 

for each depth Zo: S(Zo) = £<p (I (cp, Zo + A(Zo) cos ( (p - (po(Zo)) / N v ; with N v being the 
number of valid pixels (non-absent) along a sinusoid at depth Zo. 

10 Depth values Zo may be obtained by way of the teaching set forth in the references 

cited in the background of the invention (apparatus and methods for determining the required 
accurate depth measurements). The prior art teaches using such things as cables on a 
precision winch, a depth wheel measurement instrument, and sometimes an axial 
accelerometer whose output is integrated to help overcome the known problems inherent in 

15 precisely measuring how deep a well logging tool is in a borehole. These depth 

measurements are shown as the 100 depth unit and 110 depth unit measurements on the 
graphs in Figs. 4 and 5. 

After the averaged signals S(Zo) are calculated, and displayed as traces 23 and 24, the 
20 teaching of the prior art is used to correlate the averaged waveforms 23 and 24 to determine a 
depth offset 30 between them. One prior art technique for obtaining a measure of correlation 
to depth match two logs is disclosed in the U.S. Patent 4,312,040, issued January 19, 1982 to 
David H. Tinch, Bruce N. Carpenter and Elie S. Eliahou. The technique taught in this patent 
involves comparing two depth shifted well logs, and determining a correlation function 
25 indicating how closely the two logs match each other at a plurality of depth levels. 

The depth correlation teaching of the Tinch et al patent is succinctly summarized at 
col. 10, 1. 56 through col. 11, 1. 32. Specific detail of the correlation process taught in this 
patent is given at col. 8, 1. 16 through col. 10, 1. 39. 

30 

Briefly, referring to the summarization starting at col. 10, 1. 56 in the Tinch et al 
patent, two well logs A and B, such as from levels LI and L2, to be correlated are initially 
considered to be depth matched in accordance with the depth indicia on magnetic tape on 
which the log data is recorded, such as the depth from depth recorder 19 on recorder 16 in 
35 Fig. 1. Beginning at an initial depth level Z a in a borehole, a first assumed value of the depth 
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displacement is computed between the B log and the A log at each selected depth level from 
Z a to Z b in increments of A y . The A log is considered to be the base log and the displacement 
for the B log relative to the A log is computed. 

Concerning the depth correlation process at one depth level, the B log is effectively 
shifted one step A z at a time and a correlation function C K is computed at each such step using 
equation (1) shown in col. 4 of the Tinch et al patent. Once the correlation values of C K over 
the entire search interval 2AZc are computed, the value of K which produced the maximum 
value of Ck, designated K\ is determined. The corresponding depth displacement will be (K- 
K max /2AZ 0 . Thus, for example, if AZ D is 1 inch and (K'-K max / 2 ) is -6, the computed depth 
displacement for log B at the presently considered depth level Z will be -6 inches. However, 
this depth displacement is not computed unless C K (max) is equal to or greater than CicGim). 
This prevents weak correlation quantities from causing depth displacements. If C K (max) is 
greater than C K (Hm), the displacement SH(Z) is computed and entered into the history file for 
statistical analysis. This computed value SH(Z) is only a first assumption of the depth 
displacement which may be changed by the statistical analysis. 

This statistical analysis takes the form of declaring a depth displacement only when a 
consecutive number of identical depth displacements have been computed. By so doing, an 
occasional erroneous depth displacement caused by noise, for example, will not cause log B 
to be depth displaced. 

After a depth matching offset has been calculated for waveforms 23 and 24, according 
to the method of the above cited Tinch et al patent, the calculated offset is applied to traces 
22a - 22d shifting (depth translating) them upward as seen in Fig. 5. The shift is seen when 
comparing these traces in Fig. 5 with the same traces in Fig. 4. Traces 21a- 2 Id are thus 
depth matched to traces 22a - 22d. 

When there are more than two sets of data to be depth matched the process is 
basically the same as described above for two sets of data. The difference is that one set of 
data (from the sensors at one depth) is used as a reference set of data and the offset therefrom 
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5 to each of the other data sets is determined. The offset calculated for each of the other data 
sets is then applied to the other data sets to depth match to the reference data set. 

In Figure 6 is shown a block diagram of the steps involved in performing depth 
correlation of different sets of signal traces in accordance with the teaching of the first 
10 embodiment of the invention. At block 31 borehole images represented by the traces 21a - 
2 Id from sensors S la - Sid are obtained and oriented azimuthally. At block 32 borehole 
images represented by the traces 22a - 22d from sensors S2a - S2d are obtained and oriented 
azimuthally. 

15 As previously described, the two sets of traces 21a - 21d and 22a - 22d reflect dip of 

bedding layers in the borehole and also reflect the azimuths of the dips. At block 35 the dip 
indicated by each of these sets of images is calculated using a prior art dip calculation 
technique such as the one applying the Hough transform described in the patent mentioned in 
the background of the invention in connection with the Hough transform. 

20 

At block 36 the calculated dip for traces 21a-21dis used along with the traces 
themselves to generate a single trace 23. Similarly, the calculated dip for the traces 22a - 22d 
is used along with the traces themselves to generate a single trace 24. 

25 At block 37 a prior art single log trace depth matching technique, such as taught in the 

above cited and described Tinch et al patent, is applied to traces 23 and 24 to determine the 
optimum depth offset necessary for precisely matching the borehole traces 22a - 22d to traces 
21a - 2 Id. At block 38 the traces 22a - 22d are shifted relative to traces 21a - 2 Id to depth 
match them. 

30 

In a second embodiment of the present invention, two-dimensional well log traces 21a 
- 2 Id and 22a - 22d of Fig. 4 are transformed into three dimensional Hough space using the 
Hough transform to get three dimensional images of which only one representative slice 39 
and 40 are shown in Figure 7. The Hough transform is well known and used in the prior art 
35 to transform two-dimensional well logs into three-dimensional Hough space images from 
which dip and azimuth is determined. The Hough transform as used herein transforms a 
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5 sinusoidal waveform into a point in a three dimensional space and a plurality of sinusoidal 
waveforms making up the two-dimensional space into the three dimensional Hough space 
images. These images are made up of a plurality of image slices 39 and 40. Slice 39 is one 
of the slices of the three-dimensional image produced by the Hough transform from well log 
traces 21a - 2 Id. Slice 40 is one of the slices of the three-dimensional image produced by the 
10 Hough transform from well log traces 22a - 22d. 

Depth matching is accurately and automatically accomplished using the three- 
dimensional images created by the Hough transform from well logging traces 21a - 2 Id and 
22a - 22d. The three-dimensional Hough images have properties that make them suitable for 

15 depth matching of well logging images. Those properties are: (1) insensitivity to gaps 
(missing data) in the original well logging images (the gaps in the original well logging 
images are the spaces between traces 21a - 2 Id and between 22a - 22d in Figs. 4 and 5), (2) 
insensitivity to rotation and orientation of the pads, and (3) dip in the well logging images 
does not have to be calculated before depth matching of the well logging images is 

20 performed. 

Then, according to the second embodiment of the present invention, extremum curves 
42 and 43, representing traces 21a - 2 Id and 22a - 22d respectively, are produced from the 
three-dimensional Hough images. The extremum curves 42 and 43, reflect certain rock 
25 features in the earth, and represent peak values in the three dimensional space at a given 

depth. The image signal peaks (local maxima) in the three dimensional Hough images 39 and 
40 appear as peaks in the extremum curves 42 and 43. 

Automatic depth matching is performed directly on extremum curves 42 and 43 
30 without the need for calculating dip and azimuth of bed boundaries in an earth formation and 
compensating for it before depth matching may be performed, as in the prior art. A prior art 
depth matching algorithm is then applied to the two-dimensional extremum curves 42 and 43 
derived from the Hough image slices, such as 39 and 40, to thereby determine a depth shift 
offset that will correlate extremum curves 42 and 43. The depth shift offset that is 
35 determined is thereafter applied to traces 22a - 22d to depth match them to traces 21a - 2 Id 
as seen in Fig. 5. 



- 17- 



5 Prior art examples of depth matching algorithms that can be used to correlate the 

extremum curves may be the same as those utilized by the first embodiment and may be 
found in: (1) U.S. Patent 4,320,469, issued March 12, 1982 to William J. Frawley, and Philip 
A. Mongelluzzo, and (2) U.S. Patent 4,312,040, issued January 19, 1982 to David H. Tinch, 
Bruce N. Carpenter and Elie S. Eliahou. 

10 

The Hough transform utilized in the second embodiment of the invention is described 
in two papers: (1) J. Illingworth and J. Kittler, "A Survey of the Hough Transform", 
Computer Vision, Graphics, Image Processing 44 (1988), p. 87-116; and (2) V. F. Leavers, 
"Which Hough Transform?", CVGIP: Image Understanding 58 (1993), No. 2, p. 250-264. 
15 See also U.S. Patent 5,960,371 issued September 28, 1999 to Naoki Saito, Nicholas N. 
Bennett and Robert Burridge cited in the background of the invention; and a paper by J. 
Illingworth and J. Kittler, "A Survey of the Hough Transform" Computer Vision, Graphics 
and Image Processing, vol. 44, (1988) pp. 87-1 16. 

20 How the Hough Transform can transform images from a defined coordinate system to 

a different coordinate system is more fully explained in an article entitled "Use of the Hough 
Transform to Detect Lines and Curves in Pictures", Duda, R. P. and Hart, P. E., ACM, vol. 
15, no. 1 pp. 1 1-15. The Hough transform technique was initially developed to detect straight 
lines in binary images. However, due to its simplicity, the technique was extended to detect 

25 other types of simple analytic shapes like circles and ellipses and is used today to recognize 
patterns even when they do not have an analytic description. 

The use of the Hough transform to automatically extract features or geometric 
parameters from borehole images was first taught in U.S. Patent No. 3,069,654 issued Dec. 
30 18, 1962 to P.V.C. Hough. By using the Hough transform to analyze "dip" events they 
appear as image signal peaks in Hough three-dimensional space. The use of the Hough 
transform to measure dip and azimuth of geological beds is also taught in U.S. Pat. No. 
5,162,994, issued Nov. 10, 1992 to David O. Torres. 

35 The Hough transform makes it possible to determine, from an image, the specific 

parameters characterizing a geometrical shape such as a straight line, a circle, an ellipse, or a 
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5 sinusoid curve, such as curves 26 - 29 in Fig. 4. The Hough transform then projects points of 
these shapes into three dimensional parameter space referred to as Hough space. More 
specifically, when the Hough Transform is used to process data points on a two-dimensional 
line or curve they are, mathematically speaking, transformed into three dimensional 
parameter space where the dip and strike angles are dimensions within the parameter space. 

10 

Generally, according to Hough formalism every sinusoid is a two-dimensional curve 
that may be mathematically expressed as follows. An abscissa (p, and an ordinate Z where Z 
can be expressed as: Z = Zo + A cos ((p - 90 )• When mathematically developing this 
equation for the Hough transform we get: Z = Zo + A cos (p 0 cos (p + A sin (p 0 sin (p. We 

15 then set X 0 = +A cos (po and Y 0 = +A sin (p 0 and the equation becomes Z=Zo + Xocos(p + Y 0 
simp. In these equations A is the amplitude of a dip indicated by the two dimensional 
sinusoid in a well log image, (po is the offset (phase) of azimuth of the dip, and Zo is the zero 
crossing depth of the sinusoid. The Greek symbols cp and Z are the coordinate system of the 
logging image space, and X 0 , Y 0 and Zo are the coordinate system of the three dimensional 

20 Hough image space. 

For example, by performing a summation along the sinusoidal waveform of the 
logging image pixel measurements, a value may be obtained and associated to a point (X 0 , 
Yo, Zo) in the Hough space. The value obtained is high if the sinusoid associated to a point in 
25 the Hough space lies on a sinusoidal feature having a high pixel value in a well logging 

image. Conversely, the obtained value is low if the sinusoid associated with a point lies on a 
sinusoidal feature having a low pixel value in a well logging image. In the case where there 
is no clear sinusoidal feature in the well logging image, a medium background measurement 
is obtained. 

30 

The following algorithm is then used to generate extermum curves from three 
dimensional Hough images in order to apply a prior art depth matching algorithm: 

Let H(X 0 , Y 0 , Zo) be a 3D array of the Hough transform. 
35 1. Compute H aV erage, the mean value of H(Xo, Y 0 , Zo) one time for the whole array 
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5 2. Centralize H(X 0 , Y 0 , Zo) in the value domain: 

H c (Xo, Yo, Zo) = H(Xo, Yo, Zo) - H ave rage 

3. For each depth Zo, find the maximum absolute value V(Zo) of H c : 
V(Zo) = maxxavo { abs(H c (X 0 , Y 0 , Zo)) at Zo } 

10 From the property of the Hough transform described above one can see that, at each 

depth Zo, where the curve V(Zo) is at a local maximum, there must be a dip event (either peak 
or trough) in the form of a sinusoid, such as sinusoids 26 - 29 in Figure 4. Thus, the 
extremum curve V(Zo) characterizes the logging image and is a good candidate for applying a 
depth matching algorithm. 

15 

In Figure 8 is shown a block diagram of the steps involved in translating signals taken 
from the two sets of sensors S la - Sid and S2a - S2d into three dimensional Hough space 
(Fig. 7) to determine a depth matching offset to be used to match traces 21a - 22d and 22a - 
22d in accordance with the teaching of the second embodiment of the invention. At block 44 
20 traces 21a-21din level LI are obtained, and at block 45 traces 22a - 22d in level L2 are 

obtained, both using prior art apparatus and techniques as described in the background of the 
invention. 

At block 46 traces 21a-21din Level LI are processed using the Hough transform to 
25 produce the three-dimensional Hough image of which one slice is shown as image 39 in 
Figure 7. Also, at block 47 traces 22a - 22d in Level L2 are processed using the Hough 
transform to produce the three-dimensional Hough image of which one slice is shown as 
image 40 in Figure 7. 

30 At block 48 the three-dimensional Hough image produced in block 46 is used to 

compute a two-dimensional extremum curve, and at block 49 the three-dimensional Hough 
image produced in block 49 is used to compute a two-dimensional extremum curve. 

At block 50 a depth shift offset is calculated using the two extremum curves. As 
35 previously described, extremum curves 42 and 43 are created from the three-dimensional 

Hough images 39 and 40. An automatic depth matching calculation is performed directly on 



-20- 



5 the extremum curves, without the need for calculating dip and azimuth of bed boundaries in 
an earth formation, to calculate a depth match offset. 

At block 5 1 the calculated depth match offset calculated from the extremum curves is 
then used to depth match traces 21a - 2 Id and 22a - 22d as shown in Fig. 5. Being as signal 
10 traces 21a - 21d and signals traces 22a - 22d are offset by forty-five degrees as previously 
described, the overlaid, depth matched traces now show combined data every forty-five 
degrees around the borehole as represented in Fig. 1A. 

When there are more than two sets of data to be depth matched the process is 
1 5 basically the same as described for two sets of data. The difference is that one set of data 
(from the sensors at one depth) is used as a reference set of data and. The vertical offset 
therefrom to each of the other plurality of data sets is determined and then applied to the 
other data sets to depth match to the reference data set. As applied to the use of the Hough 
transform, a plurality of two-dimensional well log traces are each transformed into three- 
20 dimensional Hough space using the Hough transform to get three dimensional images. 
Extremum curves are then derived from each of the three-dimensional Hough images. 
Automatic depth matching is then performed directly on the extremum curves, as previously 
described. Using the extremum curves resulting from one set of data as a reference, a prior 
art depth matching algorithm is then used to compare the extremum curves representing the 
25 reference data set with the extremum curves representing each of the plurality of other data 
sets to thereby determine a depth shift offset that will correlate them to the reference set. The 
depth shift offset that is determined using each set of extremum curves is then applied to the 
corresponding data set to depth match it to the reference data set. 

30 The above description describes two embodiments of a system wherein borehole 

logging data is collected by sensors on a logging tool passing through a borehole. The sensor 
signals representing the data are transmitted via a cable to the earth surface, as shown in Fig. 
1, to be stored and processed. However, in an alternative embodiment of the invention, 
computer processor and storage capability may be added to the logging tool and depth match 

35 processing of the data from the sensors is done in real time in-situ. The depth matched data is 
then either extracted from storage after the logging tool is extracted from the borehole, or the 
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5 real time depth matched data is transmitted to the surface of the earth over cables that 

transmit data from sensors to above ground equipment as shown in Fig. 1. Alternatively, the 
depth matched data from the processor at the logging tool may be transferred from the 
logging tool to the surface using remote telemetry capability. A transmitter with the logging 
tool transmits the depth matched data determined by the receiver downhole to a receiver on 
10 the surface of the earth without the need for a physical connection between the downhole 

logging tool and the receiver. Once the depth matched data is received by the receiver it may 
be transmitted to a location from the borehole for analysis or other processing. 

What has been described above also applies to images of borehole walls and borehole 
15 core, irrespective of the measured parameter which forms the images such as, for example, 
rock conductivity, resistivity, magnetic susceptibility, photoelectric factor, density of the 
formation, amplitude of acoustic reflections, photographic brightness, etc. In addition, an 
image from one logging run in a borehole may be compared to an image from another 
logging run in the same borehole, and an image from a borehole may be compared to an 
20 image of a core from the borehole. In addition, images obtained from a cased borehole may 
be depth matched using the teaching of the present invention. Further, the present invention 
may be used to depth match logs obtained by a logging tool from a borehole during coring as 
taught in U.S. Patent 6,003,620 issued December 21, 1999 to Mukul M. Sharma, Roger T. 
Bonnecaze and Bernard Zemel. 

25 

For the sake of simplicity only, a simple arrangement for a logging tool is described. 
There are pads and sensors on two levels, each level having four pads spaced ninety degrees 
apart, and there being one sensor on each pad. The pads in one level are rotated forty five 
degrees with respect to the pads in the other level. There may be more than four pads in each 
30 level; the pads may be spaced other than equiangularly; the pads in one level may be rotated 
at angles other than forty-five degrees with respect to the pads in the other level; and there 
may be many more transducers on each pad. In addition, data or images from any type of 
transducer, sensor, scanner or camera may be processed and depth matched using the 
teaching of the present invention in both of its embodiments. 
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